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A Leading Log Approximation for Inflationary Quantum Field Theory 

R. P. Woodard a * 

a Department of Physics, University of Florida, Gainesville, FL 32611, United States 

During inflation explicit perturbative computations of quantum field theories which contain massless, non- 
conformal fields exhibit secular effects that grow as powers of the logarithm of the inflationary scale factor. 
Starobinskii's technique of stochastic inflation not only reproduces the leading infrared logarithms at each order 
in perturbation theory, it can sometimes be summed to reveal what happens when inflation has proceeded so long 
that the large logarithms overwhelm even very small coupling constants. It is thus a cosmological analogue of 
what the renormalization group does for the ultraviolet logarithms of quantum field theory, and generalizing this 
technique to quantum gravity is a problem of great importance. There are two significant differences between 
gravity and the scalar models for which stochastic formulations have so far been given: derivative interactions 
and the presence of constrained fields. We use explicit perturbative computations in two simple scalar models to 
infer a set of rules for stochastically formulating theories with these features. 



1. Overview 

This conference report represents work done in 
collaboration with Nikolaos Tsamis of the Univer- 
sity of Crete. In the next section I review the two 
conditions (masslessness and lack of conformal in- 
variance) for quantum effects to be enhanced dur- 
ing inflation. In section 3 I introduce the notions 
of infrared logarithms and the leading logarithm 
approximation. Section 4 shows that the field op- 
erator behaves like a classical random variable 
in the leading logarithm approximation. I also 
derive the approximate field equation obtained 
by Starobinskii for a massless, minimally cou- 
pled scalar with non-derivative self-interactions. 
Sections 5 and 6 generalize Starobinskii's tech- 
nique to theories which possess derivative inter- 
actions and constraints, respectively. Our results 
are summarized in section 7. 

2. Quantum Effects during Inflation 

A locally de Sitter geometry provides the sim- 
plest paradigm for inflation. To see why, con- 
sider a general homogeneous, isotropic and spa- 
tially flat geometry, 



ds* 



-dt 1 



a 2 (t)dx ■ dx 



(1) 



Derivatives of the scale factor a(t) give the Hubble 
parameter H(t) and the deceleration parameter 



H(t) 



(2) 
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The nonzero components of the Riemann tensor 
are, 



if 



= H 



(3) 



Inflation is defined as positive expansion (H(t) > 
0) with negative deceleration (q(t) < 0). On the 
other hand, stability — in the form of the weak 
energy condition — implies q(t) > — 1. At the 
limit of q = — 1 we see from © that the Rie- 
mann tensor assumes the locally de Sitter form, 

R P a»v = iPfagav-Kfrl*) ■ lt foll0WS fro m © 

that the Hubble parameter is actually constant, 
and that the zero of time can be chosen to make 
the scale factor take the simple exponential form 
we shall henceforth assume, 



de Sitter inflation 



a(t) 



(4) 
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The homogeneity of spacetime expansion in Q 
does not change the fact that particles have con- 
stant wave vectors k, but it does alter what these 
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means physically. In particular the energy of a 
particle with mass m and wave number k = \\k\\ 
becomes time dependent, 



E(t, k) 



(k/a(t) 



(5) 



This results in an interesting change in the ener- 
gy-time uncertainty principle which restricts how 
long a virtual pair of such particles with ±fc can 
exist. If the pair was created at time t it can last 
a time At given by the integral, 



t+At 

dt'E{t',k) 
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Just as in flat space, particles with the smallest 
masses persist longest. For the fully massless case 
the integral is simple to evaluate, 



t+At 

dt'E{t',k) 



1 



— HAt 



Ha(t) 



(7) 



We therefore conclude that any massless virtual 
particle which happens to emerge from the vac- 
uum with k < Ha(t) can persist forever! 

Most massless particles possess conformal in- 
variance. The change of variables dr] = dt/a(t) 
defines a conformal time in terms of which the 
invariant element is just a conformal factor 
times that of flat space, 



ds 2 



-dt 2 +a 2 (t)dx-dx = a 2 (-drj 2 +dx-dx*j .(8) 



In the (?7, x) coordinates conformally invariant 
theories are locally identical to their flat space 
cousins. The rate at which virtual particles 
emerge from the vacuum per unit conformal time 
must be the same constant — call it T — as in 
flat space. Hence the rate of emergence per unit 
physical time is, 
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It follows that, although any sufficiently long 
wavelength, massless and conformally invariant 
particle which emerges from the vacuum can per- 
sist forever during inflation, very few will emerge. 

Two kinds of massless particles do not possess 
conformal invariance: minimally coupled scalars 



and gravitons. To see that the production of these 
particles is not suppressed during inflation note 
that each polarization and wave number behaves 
like a harmonic oscillator with time dependent 
mass and frequency, 
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i(t) 
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The Heisenberg equation of motion can be exactly 
solved, 



k 2 

q + 3Hq + — 



q = 



q(t) = u(t,k)a + u*(t,k)a' < , (11) 



where the mode functions and commutation rela- 
tions are, 



u(t, k) 



and 




exp 



ik 



Ha{t) 



(12) 



The (co-moving) energy operator for this system 
is, 

E{t) = \m{t)q 2 {t) + lm(t)u; 2 (t)q 2 (t) . (13) 

Owing to the time dependent mass and frequency, 
there are no stationary states for this system. At 
any given time the minimum eigenstate of E{t) 
has energy ^oj(t), but which state this is changes 
for each value of time. The state which is 
annihilated by a has minimum energy in the dis- 
tant past. The expectation value of the energy 
operator in this state is, 



E(t) 







1 a 3 {t)\u(t,k)\ 2 + la{t)k 2 \u(t,k)\ 
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If we think of each particle having energy k/a(t) it 
follows that the number of particles with any po- 
larization and wave number k grows as the square 
of the inflationary scale factor, 



N(t,k) 



Ha(t) 
2k 



(15) 
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Quantum field theoretic effects are driven by 
essentially classical physics operating in response 
to the source of virtual particles implied by quan- 
tization. On the basis of i|15[l one might ex- 
pect inflation to dramatically enhance quantum 
effects from massless, minimally coupled scalars 
and gravitons, and explicit studies over a quar- 
ter century have confirmed this. The oldest re- 
sults are of course the cosmological perturbations 
induced by scalar inflatons 1 and by gravitons 
|5] . More recently it was shown that the one loop 
vacuum polarization induced by a charged, mass- 
less and minimally coupled scalar causes super- 
horizon photons to behave like massive particles 
in some ways |3I4I5| . Another recent result is 
that the one loop fermion self-energy induced by a 
massless and minimally coupled Yukawa scalar re- 
flects so much particle production that the end of 
inflation should reveal a degenerate gas of super- 
horizon fcrmions 

3. Infrared Logarithms and the Leading 
Log Approximation 

Two higher loop results are of special inter- 
est to us. The first is that the gravitational 
back-reaction from the inflationary production of 
gravitons induces an ever greater slowing in the 
expansion rate 7 8 9. The Lagrangian is just 
that of general relativity with a positive cosmo- 
logical constant, 



1 



16ttG 



i?-2A)V^ with A = 3H 2 .(16) 



What was actually computed |S] is the graviton 
one-point function, about a locally de Sitter back- 
ground, in the presence of a state which is free 
Bunch-Davies vacuum at t = 0. However, if the 
resulting distortion of the background geometry 
was viewed in terms of an effective energy density 
and pressure the results would be, 



P(t) = 



Pit) 



A {nH) 2 H 4 



8ttG 



2 6 7T 4 



xj-iln 2 ( a ) + o(ln(a)) 
(kH) 2 H 4 



+ 0{n% (17) 



A 

~8^G 



x| + iln 2 (a) + 0(ln( a ))| + 0(k 4 ). (18) 
(Here k 2 = 167rG.) Of course the tree order 



-) term is just the cosmological constant. 



The one loop (~ n ) effect from the kinetic ener- 
gies of inflationary gravitons is a constant |ll)lllj 
that must be subsumed into A in order for the 
universe to begin inflating with Hubble parame- 
ter H. The next order effect is secular because 
the kinetic energies interact with the total gravi- 
ton field strength which grows as more and more 
gravitons are produced. The induced energy den- 
sity of this term is negative because the interac- 
tion is attractive. 

The second result is that the inflationary parti- 
cle production of self-interacting scalars induces a 
violation of the weak energy condition on cosmo- 
logical scales |12I13| . The Lagrangian is that of a 
massless, minimally coupled scalar with a quar- 
tic self-interaction in a non-dynamical, locally de 
Sitter background, 



£ = --jpn<pdu<pg 1 ' 



- A 4 



-fJ 



(19) 



This theory can be renormalized so that, when 
released in free Bunch-Davies vacuum at t = 0, 
the energy density and pressure are |12I13| . 



P(t) 



XH 1 



A 
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+A£ln(a) + O 



0(X 2 ), (20) 
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2 6 7T 4 



3+Ae] ln(a) + o(o _1 ) [ + 0(X 2 ). (21) 



(The parameter A£ is related to the finite part of 
the conformal countertcrm.) The one loop effect 
from the kinetic energies of inflationary scalars 
is again a constant that must be subsumed into 
A. The next order effect is secular because the 
ip 4 self-interaction involves the total scalar field 
strength which grows as more and more scalars 
are produced. The induced energy density is pos- 
itive, for A > 0, because the ip 4 term adds to the 
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energy density. In this model we might also guess 
that the effect is self-limiting because the classi- 
cal restoring force tends to push the scalar back 
to zero. 

Both of these models exhibit infrared loga- 
rithms — factors of ln(a) = Ht. As inflation 
proceeds these infrared logarithms grow without 
bound until they eventually overcome the small 
coupling constants, {nH) 2 for gravity and A for 
the scalar model. One cannot conclude from 
this that there is actually a significant change 
in the background expansion rate because we do 
not know the higher order results. The legiti- 
mate conclusions are rather (1) that the expan- 
sion rate slows in the first model, and increases in 
the second, and (2) that both effects eventually 
become nonperturbatively strong. To learn what 
the outcome is at late times, after perturbation 
theory has broken down, requires a nonperturba- 
tive technique. 

One possibility for such a technique is to sum 
the leading infrared logarithms. It has been con- 
jectured that the powers of a -1 = e~ Ht — which 
are separately conserved — can be subsumed into 
renormalizations of the initial state If this 

is so then we expect the perturbative expansions 
for the energy density in both models to take the 
general form, 



P(t) 



A 



Vf( ; ' ^(coupling) 1 

2i~Z 

x ^2 G«fe( ln ( a )) 



2i-2-k 



(22) 
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where the Cik 's are pure numbers, and "coupling" 
is (kH) 2 for gravity and A for the scalar model. 
The leading logarithm approximation would be, 



p(t) loading 
logarithm 



-H 



A 

oo 

4 Cm ( coupling In 2 (a) 



(23) 



Of course other quantities besides the energy 
density exhibit infrared logarithms. For example, 
they have been evaluated at the one an d two 
loop 14 orders in the self-mass-squared of the 
scalar model l|19fl ■ Because the leading logarithm 



approximation seems to be a well-defined concept 
for any Green's function we might suspect that it 
can be implemented at the level of the Heisenberg 
field equations. 

Starobinskii has long argued that his stochas- 
tic field equations ^2] should recover the lead- 
ing infrared logarithms of any Green's function 
at each order in perturbation theory. In 1994 he 
and Yokoyama JHI were able to obtain an explicit 
late time solution of these stochastic equations for 
(|19l) which agrees with the intuitive expectation 
that the growth of its field strength is eventually 
halted by the classical restoring force of the po- 
tential. The purpose of this paper is to generalize 
Starobinskh's technique so that it can be applied 
to models with derivative interactions and con- 
strained fields such as gravity. 

4. The Physics of Infrared Logarithms 

The origin of the infrared logarithms in expres- 
sions (|17I18|) and (|2UI21I) can be understood phys- 
ically as well as from the mathematics of per- 
turbation theory. Although the physical under- 
standing is vastly more important in guiding the 
generalization we must make, it is crucial to see 
that the formalism provides an invariant separa- 
tion between infrared and ultraviolet degrees of 
freedom. We accordingly begin by explaining how 
perturbative computations are done invariantly in 
a locally de Sitter background. 

We employ dimensional regularization in posi- 
tion space. Of course the vertices are straight- 
forward to obtain in any background, so the only 
difficulty comes in finding the propagators in ar- 
bitrary spacetimc dimension D. These propaga- 
tors are expressed in terms of a de Sitter invariant 
length function we call y(x; x'), 



y(x; x') = aa' 



H- 



1 \ 2 

a 1 



.(24) 



Its physical meaning in terms of the invariant 
length £(x;x') between x^ and x ,fl is, 



y(x; x') = 4 sin 2 



-Hl(x;x') 



(25) 



For example, the propagator of a massless, con- 
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formally coupled scalar is [T7] 

H D - 2 „fD .\fl 

(47T) 



iAcf (x;x') = ^r(^-l)(l)' 
7t)~ v 2 /Vz// 



(26) 



Because y(x;x') vanishes when x' M — x M , and 
because one always interprets D so that zero is 
raised to a positive power in dimensional regular- 
ization, the coincidence limit of the conformally 
coupled propagator vanishes, 



«A cf (x; x) = . 



(27) 



Of course conformal invariance suppresses in- 
teresting quantum effects but iA c f (x; x') is useful 
for expressing the propagator of the minimally 
coupled scalar |12I13| . 

• a / / \ . . , # D - 2 r(D-i) f d 

iA^(x; x ) = iA c { (x; x ) + 



r 2 (#) / 4 



r(D-l)V!/ 



(4tt)# r(f) \D-4 
7rcot^£»^ + ln(ao')> 



H 



D-2 



1 T(n+D-1) (ij 



1 r(n+f + 1) /yxn-#+2| 



n-£ + 2 I>+2) U 



(28) 



This expression might seem daunting but it is ac- 
tually simple to use in low order computations 
because the infinite sum on the final line vanishes 
in D = 4, and each term in the series goes like 
a positive power of y(x;x'). This means that the 
infinite sum can only contribute when multiplied 
by a divergence, and even then only a small num- 
ber of terms can contribute. 

The explicit factor of ln(aa') at the end of the 
second line in (|28|l is one source of infrared log- 
arithms. It gives the secular dependence of the 
coincidence limit |18ll9l2li| . 



iA^(x;x) 



H 



D-2 



(47r)~ 

X ^l^f 7rCOt GH + 2 M«)}- (29) 

This factor of ln(aa') is also interesting in that 
it breaks de Sitter invariance. Of course <p(x) 



transforms like a scalar; the breaking derives from 
the state in which the expectation value of the 
two free fields is taken. Allen and Folacci long 
ago proved that the massless, minimally coupled 
scalar fails to possess normalizable de Sitter in- 
variant states p21 ■ A final point is that the factor 
of ln(aa') is actually an invariant, it just depends 
upon the initial value surface at which the state is 
released. In fact it is the Hubble constant times 
the sum of the invariant times from x M and x' M to 
this initial value surface. The physical reason for 
such a term is the increasing field amplitude due 
to inflationary particle production, as discussed 
in the introduction. 

The graviton propagator is much more compli- 
cated on account of its tensor indices. However, 
it also involves the A-type propagator [22l23j so 
it can also induce infrared logarithms. Any in- 
frared logarithms which appear in one loop dia- 
grams involving scalars and/or gravitons can only 
have come from this ^4-type propagator. How- 
ever, another mechanism can produce infrared 
logarithms in higher loop results such as H17I18J) 
and (|2()I21(I . This other mechanism is the growth 
of the invariant volume of the past light-cone from 
the observation point back to the initial value sur- 
face, 



V PLC (t) = J dt'a' - 1 I d D -V6» -y(x;x') 



p-i)r(^i) 



ln(a) + 0{l) 



(30) 



Factors of this quantity arise naturally whenever 
undifferentiated propagators connect an interac- 
tion vertex (at x' M ) with the expectation value of 
some observable (at x M ). 

It should be noted that one must use the 
Schwinger-Keldysh formalism [23 to obtain true 
expectation values for cases, such as cosmology, 
in which the "in" and "out" vacua either do not 
agree or are not even well defined 25 . In the 
Schwinger-Keldysh formalism the only net con- 
tributions arise from interaction vertices which lie 
on or within the past light-cone of some observa- 
tion point. The more familiar in-out matrix ele- 
ments would harbor virulent infrared divergences 
from integrating over the exponentially large vol- 
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lime of the inflationary future The causal- 

ity of the Schwinger-Keldysh formalism regulates 
these infrared divergences but simple correspon- 
dence implies that the regulated expressions must 
grow without bound at late times. 

To summarize, infrared logarithms derive joint- 
ly from explicit factors of ln(aa') in the propaga- 
tors and from integrating undifferentiated prop- 
agators over the past light-cone. Our reason for 
discussing the mathematics of perturbation the- 
ory is to emphasize that the ultraviolet regulariza- 
tion is invariant. This has a crucially important 
consequence for understanding the physical origin 
of infrared logarithms: it means that there should 
be a constant dynamical impact from all modes 
whose physical wavelength ranges from zero (the 
far ultraviolet) to any fixed value. 

Recall from section 2 that quanta are labeled 
by constant wave vectors fc, and that the mode 
with wavenumber k = \\k\\ begins to experience 
significant inflationary particle production when 
k/a(t) = H. This suggests that we distinguish 
"infrared" and "ultraviolet" modes on the follow- 
ing basis, 

Infrared H < k < H a(t) , (31) 

Ultraviolet =*> k > H a(t) . (32) 

The ultraviolet modes are those whose physical 
wavelengths (= 2ita(i)/k) range from zero to the 
constant 2tt/H. With any invariant regulariza- 
tion the dynamical impact of such modes must 
be constant because they lie in an invariantly de- 
fined range. 

Had we taken the lower limit in expression Il31|l 
to k = the infrared phase space would also have 
extended over a constant physical range. How- 
ever, taking k down to zero has long been known 
to result in infrared divergences We regu- 

late these by working on T ^ 1 x R, with the co- 
ordinate toroidal radii equal to 2ir/H [2EJ. (We 
also typically make the integral approximation to 
discrete mode sums.) There are other techniques 
|28| but they all cause the effective infrared phase 
space to increase as the universe inflates. This 
growth is the physical source of infrared loga- 
rithms. Essentially what happens is that the av- 
erage field strength increases as it receives contri- 
butions from more and more infrared modes. If 



any interactions involve the undifferentiated field 
then this growth can enter physical quantities. 

These assertions are easy to confirm in the con- 
text of the scalar ip 4 model ()19f) in D = 3 + 1 
dimensions. The field equation is, 



V 2 A 

Cp + ZHip -<p+-tp 3 = 

a z 6 



(33) 



The full (ultraviolet plus infrared) perturbative 
initial value solution is obtained by iterating the 
following equation, 

A f* 

<p(t,x) =ip (t,x)-- / dt'a 3 (t') 

x J d 3 x'G iet (t,x;t' ,x*^y 3 {t' ,x'). (34) 
The retarded Green's function is, 

Gret (t) X] if) x'j 




• (35) 

v IX U ' ✓ I 

The free field is, 

(*.*) = J^{e llS u(t,k)a(k)+c.c.}, (36) 



where the mode function and commutation rela- 
tions are, 



u(t, k) 



H 



1 



ik 



and 



a(k),a\k')] = (2nf6 3 (k-k / y (37) 

Note that although ipo(t,x) is only the order A 
part of the solution, it and its first derivative 
agree exactly with the full field on the initial value 
surface. 

To excise the ultraviolet modes l|32|l we iterate 
what is essentially the same equation, 

A 

$(t,x) = $ (t,x)- ~ J dt'a 3 (t') 

x I d 3 x'G Tet (t,x-,t',x')<5> 3 (t',x') , (38) 
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but with the zeroth order solution restricted to 
only infrared modes, 



(27T) 

^ S u{t, k)a{k) 



(t,k)a\k)}, (39) 



This model is completely free of ultraviolet diver- 
gences, so we are justified in taking D = 3 + 1. 
To see that the model also reproduces the lead- 
ing infrared logarithms at tree order it suffices to 
take the expectation value of $§, 



$g(t,£) 



H 2 c 



,(40) 



(41) 



Comparison with (|29() (for D = 3 + 1) reveals 
exact agreement between the ln(a) terms. 

Our purely infrared field operator <I>(t, x) does 
not quite obey the field equation (|33(l because the 
kinetic operator fails to annihilate &o(t, %)■ When 
time derivatives act upon the time dependent up- 
per limit of the mode sum, 

V 2 

= T{t,x')+G(t,x)+?>HF{t,x), (42) 

it produces momentum space surface terms in- 
volving the sources, 

F(t,2) = Hj-^k8(k-Ha{tj) 

xje^uCt, k)a(k) + e-^ s u*(t, fc)a f (fc)}, (43) 

/d 3 k / \ 
—-kS[k-Ha(t)j 

*(t,k)a\k)}. (44) 



^—ik-x • * / 

-e u ' 



(2tt) 3 

x{e iSs u(t,k)a(k) 

Note that the mode functions and their deriva- 
tives are simply constants at k = Ha(t), 

H 

k=Ha{t) 



u(t, k) 
u(t, k) 



'2k 3 

H 2 



l-i)e\ 



k=Ha{t) 



(45) 
(46) 



In view of l|42|l the equation obeyed by Q(t, x) 
is not (|33|l but rather, 

$-^-^)+3ij($-J')-^-$ + ^$ 3 = .(47) 

This equation does not leave the original stress- 
energy tensor conserved because stress-energy is 
being continually dumped into the truncated sys- 
tem by ultraviolet modes which redshift past the 
horizon. We can account for this by modifying 
what we call T^ IV . To motivate the modification 
it is useful to write down the form we expect for 
the divergence of the stress-energy. For the zero 
component we should get, 



V = -T a 



H 



3T n , 



+ ^T 0lJ , (48) 



V\ A 3 



The spatial components should read, 

1 



a 



(49) 



(50) 



-ft* 



V\ A 3 
a z 6 



(51) 



We can enforce <|49|1 and (|51|l with a stress- 
energy of the form, 



1 
2a~ 2 



24 



-$ 4 



T 0l - ($-.F)&$, 
'/,, - di$dj$ - gij 

<t" ; +d i S j + d j S i 



2a- 



(52) 
(53) 



A 



24 



Once this form is assumed the nonlocal source 
terms of the purely spatial components are deter- 
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mined by conservation to be, 



S- — 



gdi$ + ($-F)diF 
a 4 did k r 



S L = 

s = - 



a 2 d k 



V 2 



ga fe $ + ($-^)a fc jf 



Aa 2 
6# 



(55) 
(56) 
J="$ 3 (57) 



We have therefore achieved a completely consis- 
tent model of just the infrared modes that repro- 
duces the leading infrared logarithms. 

Although the zeroth order field §>o(t,x) con- 
tains only infrared modes it is still quantum me- 
chanical in that the field and its first time deriva- 
tive do not commute, 



(27T) 

u(t, k)u*(t, k) ~ u*(t, k)u(t, k) 



4-K 2 a 3 



Ha 



dk k 



H 



2 sin(fcAx) 
kAx : 



k) 

4fe - Ai? ,(58) 
(59) 



Because it is a superposition of creation and an- 
nihilation operators <po (t, x) is still a quantum op- 
erator and hence liable to take any value. Such a 
random but commuting variable might be termed 
"stochastic." 

How should one include corrections to 4>o(t, x)? 
We could simply replace $o by <fio in Iter- 
ating such a relation would capture the leading 
infrared logarithms however, we may as well sim- 
plify the Green's function along the same lines as 
the mode function. The full Green's function can 
be expressed in terms of the full mode function 
as follows, 

G rct (t,x;t',x') =i9{t-t') 1 



(2tt) 3 



u(t, k)u*(t', k) — u*(t, k)u(t , k) 



ik-Ax 



(63) 



The long wavelength expansion of the bracketed 
term is, 



u(t,k)u*(t',k)-u*{t,k)u(t',k) 



1 



1 



(64) 



4TT 2 (aAx) 3 
sm(aH Ax) — aHAx cos(aH Ax) 



< kn{HAx) + HAxcos(HAx) ] ' ^ 

(In these formulae we define Ax = x — x' and 
Aa; = ||Ax||.) However, note that the leading in- 
frared logarithm in (|41Jl derives entirely from the 
constant first term of the long wavelength expan- 
sion of the mode function. 



H 



1 / k \ 2 i f k \ 3 
2p1 + 2 \~Ha) + 3 \~Ha) ^ 



(61) 



We would get precisely the same infrared loga- 
rithm from a free field with only this term, 



<j>o{t,x) 



d 3 k 
(2^)3 

( „ik-S 



e[Ha{t)-k 



a(k) + e" l£ -V(fc)}. (62) 



However, the field <po(t,x) is "classical" in the 
sense that it commutes with its time derivative. 



Retaining only this leading term in l|63[) gives, 

Giet (t, X] t' , x\ 



1 

3H 



9{t-t')S 3 (x- x')\ 



1 



1 



(65) 



It is actually superfluous even keeping the 1/a 3 
term, which can never result in a leading infrared 
logarithm. We will therefore recover the lead- 
ing infrared logarithms by iterating the simplified 
equation, 



<j)(t,x) = 4>o{t,x) 



A 



18H 



dt'<t) 3 {t',x) 



(66) 



It remains to infer the local differential equa- 
tion which the full stochastic field <p(t, x) obeys. 
Note that the time derivative of the stochastic 
free field is a momentum space surface term, 



4>o(t,x) = f(t,x 
r d 3 k 
W) 



(67) 

rf ^(fc-ffa)-^={e if «(fc)+c.c}, (68) 
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Hence the full stochastic field obeys, 
A 



3H 



(j>(t,sf)-f(t,Zj) + -<f> 3 (t,x) =0. (69) 



5. Derivative Interactions 

Of course ljfi9*|l is a Langevin equation with the 
white noise provided by f(t,x). For the tech- 
niques by which the associated Fokker-Planck 
equation for the probability density can be solved 
at asymptotically late times we recommend the 
excellent paper by Starobinskii and Yokoyama 
[*R)| . Our main purpose is to derive the analogous 
Langevin equation for quantum gravity. This en- 
tails generalizing the derivation to cover deriva- 
tive interactions (this section) and constrained 
fields (section 6). 

The key property of the Langevin equation we 
seek is that it should capture the leading infrared 
logarithms at each order in perturbation theory. 
Because we shall have to make guesses about the 
equation it is desirable to be able to check these 
guesses against explicit perturbative computa- 
tions in the full theory. The time required to do 
this in quantum gravity is prohibitive. We there- 
fore studied a simple scalar model with derivative 
interactions, 

C = -\d^d v Ag^^—g- ±d„Bd v Bg^V=9 



A 



A 2 d fl Bd u Bg^^. (70) 



Both fields A(t,x) and B(t,x) are massless, min- 
imally coupled scalars at zeroth order, so both 
have the A-type propagator (|28[l. With a little 
work we can evaluate the one loop corrections to 
the square of each field |2H], 



A 2 {t,x) 
B 2 (t,x) 



leading 
log 



leading 
log 



4tt 2 
4tt 2 



ln(a) + 0(A 2 ), (71) 
ln(a) 



XH 4 



ln 2 (a)+0(A 2 ). (72) 



The full field equations for this system are, 



1 5S A , 



(73) 



7=g% = B '^ + XAA ^ + ^A 2 B^ = 0.(74) 

We can solve them perturbatively by iterating the 
following system, 

A(t,x) = A fl (t,x) - ~ J dt'a' 3 J d 3 x'G xet (x;x') 

xA(x')B^(x')B>»(x'), (75) 
B(t, x) = B b (t, x)+\J dt'a' 3 J d 3 x'G lct (x; x') 

x{a(.t')4m(^)^(^) + ^ 2 (^)^(^)}- (76) 

(In these formulae we have sometimes compressed 
space and time arguments to a single 4-vector: 
(t, x) = x.) The two free fields are independent 
copies of the scalar free field 



Ut,?) = J^{e^ S u(t,k)a(k)+c.c.}, (77) 
'fr(t,2) = J-0^{e rk Mt,k)/3(k) + c.c.}. (78) 



Here a(k) and (3(k) are independent, canonically 
normalized annihilation operators and u(t, k) is 
the mode function (|37|l of the massless, minimally 
coupled scalar. 

What is the simplified stochastic analogue that 
recovers the leading infrared logarithms? From 
the analysis of the previous section it is appar- 
ent that differentiating a propagator precludes it 
from contributing an infrared logarithm. Said 
another way, it is always preferable to have a 
derivative act on the rapidly varying scale factor 
rather than the slowly varying field. We therefore 
expect that the Hubble friction term dominates 
the scalar d'Alembertian, 



. V 2 

-A-3HA+ —A 



3HA. 



(79) 



No term is comparably important when the two 
derivatives act on different fields, 



A uB'V = -AB + X;VA-VB — ► 0, 



B a B>>* = -B 2 + -^WB VB 



1 



0. 



(80) 
(81) 



The rest of the derivation is identical to that of 
the previous section. We again excise sub-horizon 
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modes and take the long wavelength limit of the 
mode function, 



A (t,x 



B (t,x 







■") 




H 






(82) 


•>-! 




-*) 




H 

V2k? 






(83) 



We also take the long wavelength limit (|(j5|> of the 
Green's function. Putting everything together 
gives the following simplified iterative equations 
for reproducing the leading infrared logarithms, 



A IR (t,x) = A (t,x), 
B 1R (t,x) = B (t,x) 
A 



dt'A 2 R (t',x)B 1R (t',x) 



(84) 
(85) 



We differentiate to obtain local equations, 

-3#(i m -/ A ) = 0, (86) 



3XH 



-A m B 1R 



-3H[B m -f B )- . } 



Here the stochastic sources are, 
d 3 k 



0. 



(87) 



(2tt) 3 

J Ak-S 



d(k-Ha(t) 



{e^a^+e-^Vffe)}, (88) 

/ <SH*-""C>) 

2fc ={e i ^(fc)+e- i ^ t (fe)}. (89) 



2fc 

2n 
H 2 



fu(t,x) 



Note the curious fact that even interactions which 
contain derivatives are free of stochastic source 
terms. 

Although we would like to check it on other 
quantities and at higher orders, equations Q86I - 
I87|) do reproduce the known leading infrared log- 
arithms I|71I72|1 . An amusing feature of this sys- 
tem is that we can obtain an explicit operator 
solution rather than just a solution for the prob- 
ability distribution at asymptotically late times! 



Of course we can read off A IR from 1)84)1 . A few 
simple rearrangements gives a closed form solu- 
tion for B, R as well, 



(90) 



It exhibits a sort of spacetime-dependent field 
strength renormalization whereby each contin- 
gent of /3-modes which experiences horizon cross- 
ing is attenuated by the factor 1 + ^A^{t,x). 

6. Constrained Fields 

Gravity also possesses constrained fields. We 
might anticipate that these possess no stochastic 
source terms because they have no independent 
degrees of freedom. Finding a scalar analogue 
model in which to check this is difficult because 
pure scalar models lack gauge constraints. How- 
ever, it is possible to go the other way by formu- 
lating gravity in a covariant gauge (3U] for which, 
in D = 4 spacetime dimensions, all the free fields 
are either massless, minimally coupled scalars or 
else massless, conformally coupled scalars. We 
can then compare scalar models involving the two 
sorts of fields. 

A simple model of the desired type is, 



C = -~d„Ad u Ag< t \ 



-9- 



The A field again has the A-type propagator (|28|) 
while the C field has the conformal propagator 
(|26l) . Taking the VEV's of the squares of each 
field gives, 



^■^ as -^Ha) + 0(^). (92) 

log 



n A 2 (t,x) 
n c 2 (t,x) 

The exact field equations are 



=0 + O(k 2 ). 

log 



(93) 



(94) 



-== y"t = A^ — kH 2 AC = 0, 
f=g 5A " 

2= ^=C>» -2H 2 C-IkH 2 A 2 =0. (95) 
f=g SC M 2 
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As usual the perturbative solution follows by in- 
terating the integrated field equations, 

A(t,x) = A b (t,x) - kH 2 J dt'a z {t') J d 3 x' 

x G rct (t, x; t', x'^j A(t', x')C{t', f),(96) 
1 



C{t,x) = C h {t,x) - -kH 2 / dt'a 3 (t') 
2 Jo 

xT Ict (t,x;t',x^A 2 (t',x'). 



(97) 



The free A-type field is identical to (|77|l . How- 
ever, the free C-type field has a different mode 
function, 

C b {t,x) = /-^{e^(t,t) 7 (t)+c.c.). (98) 



v(t,k) 



(2tt)3 



-i 1 



'2k a 



. e Ha 



(99) 



The C-type field also has a different retarded 
Green's function, 



rv 



: (t,x; t\x'^j 



H 2 
— 0(t-t 
4tt v 



aa' 'H\\x—x' 



(100) 



We can implement the leading logarithm ap- 
proximation on the ^4-type field equation using 
techniques which are by now familiar, 



A m (t,x) = A (t,x) 

1 r* 

--kH / dt'A m (t',x)C 1R {t\x). 
-J Jo 



(101) 



Achieving a similar reduction for the C equa- 
tion requires several new insights. First, the C- 
type mode function 199|) lacks the fc~ 2 singularity 
needed to produce infrared logarithms. It follows 
that the C-type free field truncates to zero in the 
leading logarithm approximation! The long wave- 
length limit of the C-type retarded Green's func- 
tion is also different, 



r r 



(t, X]t',X* ^ 



(102) 



These insights allow us to write a leading loga- 
rithm approximation for the C equation as fol- 
lows, 

Cin {t,x) = -^-fdt'(- - ^)A 2 1R {t',x). (103) 
2 J \ a a z 1 

A final simplification is to integrate by parts on 
the rapidly varying scale factors and neglect the 
slow variation in A m (t', x) to obtain a local equa- 
tion, 



Cm — ^ K A 2 R . 



(104) 



We can rewrite the simplified leading logarithm 
field equations as follows, 

- 3-ff(i ra -/ A ) - kH 2 A ir C ir = 0, (105) 
-2H 2 C m -^KH 2 A 2 lR = 0. (106) 

The C equation (|106|l could be derived from l|95J) 
by arguing that it is preferable to act deriva- 
tives on the rapidly varying scale factor rather 
than the slowly varying field. Hence the confor- 
mal d'Alembertian is dominated by the curvature 
term, 



C< 



-RC 
) 

-2H 2 C. 



-C-3HC 



-C-2H 2 C 



(107) 



We would again like to check more expectation 
values and higher orders, but our putative equa- 
tions (|105I106(I do reproduce the known leading 
infrared logarithms (|92I93|) . Substituting l|104|l 
into the A equation (|105|) gives, 



.4, 



/a 



k 2 H 
~Y2 



-At 



(108) 



It follows that A IR (t, x) and <fi(t, x) obey the same 
equation with A = ^k 2 H 2 \ 

7. Conclusions 

We can at length state what seem to be the gen- 
eral rules for inferring simplified operator equa- 
tions which reproduce the leading infrared loga- 
rithms: 
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1. At each order in the field (e.g., ip 1 , ip 2 , etc.) 
retain only the term with the smallest num- 
ber of derivatives. 

2. Each time derivative in the linear terms has 
a stochastic source subtracted. 

Of course we would get the same leading infrared 
logarithms by solving (|47|l for x). The enor- 
mous advantage of solving the stochastic equa- 
tion Q69[) instead is that the field <p(t, x] can be 
regarded as classical. Whereas there is little hope 
of being able to exactly solve the Heisenberg field 
equations for an interacting, four dimensional 
quantum field, one can sometimes solve classical 
field equations. Starobinsky and Yokoyama did 
it for a minimally coupled scalar with arbitrary 
potential [TB] . We have just seen how it can be 
done for models with derivative interactions and 
with constrained fields. The great task now is to 
see what happens for quantum gravity. 
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